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We analyze the properties of a 19 dimensional Galerkin approximation to a parallel shear flow. 
The laminar flow with a sinusoidal shape is stable for all Reynolds numbers Re. For sufficiently 
large Re additional stationary flows occur; they are all unstable. The lifetimes of finite amplitude 
perturbations show a fractal dependence on amplitude and Reynolds number. These findings are in 
accord with observations on plane Couette flow and suggest a universality of this transition szenario 
in shear flows. 

47.20.Ft, 47.20.Ky, 47.15.Fe, 05.45. +b 



Oh! 

< 

00 
(N 



> ■ 
OO 1 

^t- : 
o ■ 
■ 

o : 
o\ . 

On ■ 

: 

>v 
■ 

i 1 

o : 

^ ■ 

O ■ 



I. INTRODUCTION 

In many flows the transition to turbulence proceeds 
via a sequence of bifurcations to flows of ever increas- 
ing spatial and temporal complexity. Analytical and ex- 
perimental efforts in particular on layers of fluid heated 
from below p],[| and fluids between rotating concentric 
cylinders [^|,[5| have lead to the identification and veri- 
fication of several routes to turbulence, which typically 
involve a transition from a structureless laminar state to 
a stationary spatially modulated one and then to more 
complicated states in secondary and higher bifurcations. 

Transitions in shear flows do not seem to follow this 
pattern ||. Typically, a transition to a turbulent state 
can be induced for sufficiently large Reynolds number 
with finite amplitude perturbations, just as in a subcrit- 
ical bifurcation. However, in the most spectacular cases 
of plane Couette flow between parallel plates and Hagen- 
Poiseuille flow in a pipe j5j, there is no linear instability 
of the laminar profile for any finite Reynolds number that 
could give rise to a subcritical bifurcation. The turbulent 
state seems to be high dimensional immediately, with- 
out clear temporal or spatial patterns (unlike the rolls 
in Rayleigh-Benard flow). And the transition seems to 
depend sensitively on the initial conditions. Based on 
these characteristic features it has been argued that a 
novel kind of transition to turbulence different from the 
well-known three low-dimensional ones is at work [0 . 

Recent activity has focussed on three features of this 
transition: the non-normality of the linear eigenvalue 
problem the occurence of new stationary states 

without instability of the linear profile |l^-[l5| and the 
fractal properties of the lifetime landscape of perturba- 
tions as a function of amplitude and Reynolds number 
fbSf . The non-normality of the linear stability problem 
implies that even in the absence of exponentially growing 
eigenstates perturbations can first grow in amplitude be- 
fore decaying since the eigenvectors are not orthogonal. 
During the decay other perturbations could be amplified, 
giving rise to a noise sustained turbulence Q . The am- 
plification could also cause random fluctuations to grow 



to a size where the nonlinear terms can no longer be ne- 
glected |n],|ll| • Then the dynamics including the nonlin- 
ear terms could belong to a new asymptotic state, differ- 
ent from the laminar profile, perhaps a turbulent attrac- 
tor. Presumably, this attractor would be built around 
stationary or periodic solution. Here, the observation of 
tertiary structures ]l^-[l5[| comes in since they could form 
the basis for the turbulent state. Finally, the observation 
of fractality in the lifetime distribution suggests that the 
turbulent state is not an attractor but rather a repeller: 
Infinite lifetimes occur only along the stable manifolds 
of the repeller, all other initial conditions will eventually 
decay. Permanent turbulence would thus correspond to 
noise induced excitations onto a repeller. 

In plane Couette flow some of the features described 
above have been identified, but only with extensive nu- 
merical effort |^2| ^4[@ . The aim of the present work is 
to present a simple model that is based on the Navier- 
Stokes equation and captures the essential elements of the 
transition. It is motivated in part by the desire to obtain 
a numerically more accessible model which perhaps will 
provide as much insight into the transition as the Lorcnz 
model 1 18 1 for the case of fluids heated from below (pre- 
sumably at the price of similar shortcomings). The two 
and three degree of freedom models proposed by vari- 
ous groups (and reviewed in |is|] ) to study the effects of 
non-normality mock some features of the Navier-Stokcs 
equations considered essential by their inventors but they 
are not derived in some systematic way from the Navier 
Stokes equation. The model used here differs from the 
one proposed by Waleffe jlO) in the selection of modes. 

Attempts to built models for shear flows using Fourier 
modes immediately reveal an intrinsic difficulty: In the 
case of fluids heated from below the nonlinearity arises 
from the coupling of the temperature gradient to the flow 
field so that two wave vectors, k and 2k, suffice to ob- 
tain nonlinear couplings. In shear flows, the nonlinear- 
ity has to come from the coupling of the flow field with 
itself through the advection term (u • V)u. This im- 
poses rather strong constraints on the wave vectors. At 
least three wave vectors satisfying the triangle relation 
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ki + k2 + k3 = are required to collect a contribution 
from the advection term. A minimal model thus has at 
least six complex variables. Three of these decay mono- 
tonically to zero, leaving three for a nontrivial dynamics. 
In the subspaces investigated (B.E., unpublished), the 
most complex behaviour found is a perturbed pitchfork 
bifurcation, which may be seen as a precursor of the ob- 
served dynamics: for Reynolds numbers below a critical 
value, there is only one stable state. Above that value 
a pair of stable and unstable states is born in a saddle- 
node bifurcation. The stable state can be excited through 
perturbations of sufficient amplitude. The basins of at- 
traction of the two stable states are intermingled, but the 
boundaries are smooth. 

Thus more wave vectors are needed and they have to 
couple in a nontrivial manner to sustain permanent dy- 
namics. The specific set of modes used is discussed in 
section [nj. It is motivated by boundary conditions for 
the laminar profile and the observation that wave vectors 
pointing to the vortices of hexagons satisfy the triangle 
conditions in a most symmetrical manner. Other than 
that the selected vectors are a matter of trial and error. 
In the end we arrive at a model with 19 real amplitudes, 
two force terms and 212 quadratic couplings. Without 
driving and damping the dynamics is energy conserving, 
as would be the corresponding Euler equation (suitably 
truncated). Moreover, the perturbation amplitudes can 
be put together to give complete flow fields. Thus the 
model has a somewhat larger number of degrees of free- 
dom, but the dynamics should provide a realistic approx- 
imation to shear flows. 

The outline of the paper is as follows. In section [n] we 
present the model, in particular the selected wave vec- 
tors, the equations of motion and a discussion of symme- 
tries. In section [II we focus on the dynamical properties 
of initial perturbations as a function of amplitude and 
Reynolds number. In section |y| we discuss the station- 
ary states, their bifurcations and their stability proper- 
ties. We conclude in section [y] with a summary and a 
few final remarks. 



II. THE MODEL SHEAR FLOW 

Ideal parallel shear flows have infinite lateral extension. 
Both in experiment and theory this cannot be realized. 
We therefore follow the numerical tradition and chose 
periodic boundary conditions in the flow and neutral di- 
rection. The flow is confined by parallel walls a distance d 
apart. A convenient way to built a low dimensional model 
is to use a Galerkin approximation. Solid boundaries 
would require the vanishing of all velocity components 
and complicated Galerkin functions where all the cou- 
plings can only be calculated numerically. However, un- 
der the assumption that here as well as in many other sit- 
uations the details of the boundary conditions effect the 
results only quantitatively but not qualitatively, we can 



adopt free-free boundary conditions on the walls and use 
simple trigonometric functions as basis for the Galerkin 
expansion. Similarly, the nature of the driving (pressure, 
boundary conditions or volume force) should not be es- 
sential so that we take a volume force proportional to 
some basis function (or a linear combination thereof). 
This still leaves plenty of free parameters to be fixed be- 
low. 



A. Galerkin approximation 

We expand the velocity field in Fourier modes, 

u(x,t) = 5>(M)e ik - x . (l) 



Incompressibility demands 

u(k,t) • k = 0. 



(2) 



The Navier-Stokes equation for the amplitudes u(k, t) 
becomes 

9 t u(k, t) = - ipkk - i ^2 (u(p, t) ■ q) u(q, t) 

p+q=k 

-z/k 2 u(k,t) + / k (3) 

where pk are the Fourier components of the pressure (di- 
vided by the density) , v is the kinematic viscosity and /k 
are the Fourier components of the volume force sustain- 
ing the laminar profile. 

There are three constraints on the components u(k): 
incompressibility Q), reality of the velocity field, 



u(-k) = u(k)* 



(4) 



and the boundary conditions that the flow is limited by 
two parallel, impenetrable plates. The ensuing require- 
ment u z (x, y, z) = at z = and z = d (where d is the 
separation between plates) is most easily implemented 
through periodicity in z and the mirror symmetry 

I u y J (x, y, -z) = I u y Ux,y,z), (5) 

which in Fourier space requires 

/u x \ / u% \ 

I u y I ( k>xi ky,k z )= I Uy J (k x ,ky,k z ). (6) 

V u z ) \ -u* ) 

This is not sufficient to fix the coefficients: the dynamics 
also has to stay in the relevant subspace, and thus the 
time derivatives have to satisfy similar requirements. 
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B. The wave vectors 



n=(l,0,0) T m=(0,l,0) T . 



(11) 



The choice of wave vectors is motivated by the geome- 
try of the flow and the aim to include nonlinear couplings. 
The basic flow shall be a flow in y-direction, neutral in 
the ^-direction and sheared in the z-direction. Thus we 
take the first three wave vectors in z-direction, 



ki = 



k 2 = 



k, = 



(7) 



The negative vectors — k^ also belong to the set but will 
not be numbered explicitely. In these units, the period- 
icity in the z-direction is 27r, so that the separation be- 
tween the plates is d — tt because of the mirror symmetry 
(|). The amplitude u(ki) will carry the laminar profile 
and u(k 3 ) can be excited as a modification to the lam- 
inar profile. k 2 is needed to provide couplings through 
the nonlinear term. These three vectors satisfy a triangle 
identity ki + k 2 — k 3 = 0, but the nonlinear term vanishes 
since they are parallel. 

The next set of wave vectors contains modulations in 
the flow and neutral direction, 



k 4 = 



(0), 



k fi = 




(8) 



Together with — k^ they form a regular hexagon, so that 
they provide nontrivial couplings in the nonlinear term. 
The periodicity in flow direction is 47r/v3, in the neutral 
direction it is 4ir. 

Finally, this hexagon is lifted upwards with ki and k 2 
to form the remaining 12 vectors, 



k 7 


= ki 


+ k 4 


k 8 


= ki 


+ k 5 


k 9 


= ki 


+ k 6 


kio 


= k x 


-k 4 


kn 


= k 4 


-k 5 


ki 2 


= k 4 


-k 6 


ki3 


= k 2 


+ k 4 


ki 4 


= k 2 


+ k 5 


ki 5 


= k 2 


+ k 6 


ki 6 


= k 2 


-k 4 


ki7 


= k 2 


-k 5 


ki 8 


= k 2 


-k 6 



(9) 

The full set k,, i = 1 . . . 18 is shown in Fig. |l|. 

The Fourier amplitudes u(k^) have to be orthogonal to 
ki because of incompressibility (|2|) . If they are expanded 
in basis vectors perpendicular to k^, the pressure drops 
out of the equations and need not be calculated. We 
therefore chose normalized basis vectors 



n(ki) 
m(ki) 



h>r. k ? _ 



kyk 



U2 _|_ U2 ' V2 



(ky, k x ,0) 




k 2 J{kl 



The amplitudes of the velocity amplitude are now ex- 
panded as 



u(ki, t) = a(ki, t)n(kj) + /3(kj, i)m(kj 



(12) 



The impenetrable plates impose further constraints on 
the a(ki) and /3(kj). For i = 1, 2 and 3 the wave vector 
has no components in the x- and y-directions, so that a 
and (3 have to be real. For i = 4, 5 and 6 the veloc- 
ity field cannot have any components in the z-direction, 
hence a = 0. The remaining wave vectors k^ and — k; 
with i = 7, . . . , 18, a total of 24, divide up into six groups 
of 4 vectors each, 



k — {k x , ky , k z ) , k — ( k x , ky,k z ), k and k . 



(13) 



The groups are formed by the vectors and their negatives 
in the pairs with indices (7,10), (8,11), (9,12), (13,16), 
(14,17) and (15,18). The amplitudes of the vectors in 
the sets are related by 



a(k) 
/3(k) 



a(- 
/*(- 



-k)« 

-kr 



-a(k') 
-/3(k') 



(14) 



Thus the full model has 6 + 6 + 6x4 = 36 real am- 
plitudes. Restricting the flow by a point symmetry 
around Xo = (0, 0, rr/2) T eliminates the contributions 
from k 2 and some other components, resulting in a 19- 
dimensional subspace with nontrivial dynamics and the 
following amplitudes: 



a(ki) 
a(k 3 ) 
/?(k 4 ) 
a(k 7 ) 
a (kg) 
a (k 9 ) 
a(k 43 ) 
a(ki 4 ) 
a(k 45 ) 



2/1 
2/3 

W5 
2/8 

2/10 
2/12 

iyu 



/9(ki) = 
/?(k 3 ) = 
/3(k 5 ) = 
/9(k 7 ) = 
/?(k 8 ) = 
/?(k 9 ) = 
/3(k 13 ) 
/3(k M ) 
/3(k 15 ) 



1)2 
2/4 

: iye 

2/9 

: 2/11 

= 2/13 

= W15 

= W17 

= iyw ; 



/3(k 6 ) = iy 7 



(15) 



components not listed vanish or are related to the given 
ones by the boundary conditions (|l4|). A complete list- 
ing of the flow fields associated with the coefficients 
Ui such that u = 2/i u i as we H as °f the equations of 
motion are available from the authors. 



(10) 



C. The equations of motion 



so that n, m and k form an orthogonal set of basis vec- 
tors. For the negative vectors — ki we chose the basis 
vectors n(— k^) = n(k;) and m(— k^) = — m(ki). If the x 
and y components of k vanish, the above definitions are 
singular and replaced by 



In this 19-dimensional subspace y\ . . . y±g the equations 
of motion are of the form 



m = ^2 AjkVjyk - vK.^ + fi 



(16) 
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Of the driving force all components but and f± van- 
ish. Moreover, if the /'s are taken to be proportional to 
u, the resulting laminar profile has an amplitude inde- 
pendent of viscosity (and thus Reynolds number). These 
components give rise to a laminar profile that is a super- 
position of a cos(z) profile (from $2) and a cos(3z) profile 
(from fi). This allows us to approximate the first two 
terms of the Fourier expansion of a linear profile with 
velocity u v = ±1 at the walls, 



1 



u 



(cosz + - cos3z) e y 



(17) 



that can be obtained with a driving /a = Ais/tt 2 and 
/ 4 = 4^/97r 2 (see Fig. |). 

The nonlinear interactions in the Navier-Stokes equa- 
tion conserve the energy E — \ J dVu 2 . In the 19- 
dimensional subspace, the corresponding quadratic form 
is 



19 



E = v E^+ 2 E^ 



(18) 



The above equations conserve this form without driving 
and dissipation. With dissipation but still without driv- 
ing, the time derivative is negative definite, indicating a 
monotonic decay of energy to zero. 

Finally, we define the Reynolds number using the wall 
velocity of the linear profile, uq = 1, the half width of 
the gap, D = d/2 = ir/2 and the viscosity v, 



Re = uqD jv = 7r /2v . 



(19) 



The other geometric parameters are a period 47r/v3 in 
flow direction and Att perpendicular to it. 



D. Symmetries 



/ U x \ / -U x \ 

Ti : I u y J (x, y, z) ->• I u y \(-x,y,z). (22) 
\u~ / \ u z ) 



and two shifts by half a lattice spacing, 



T 9 



T, 




(x + 2ir,y,z) 



(23) 



(x + 7T, y + tt/V3, z) . (24) 



When applied to the flow these transformations induce 
changes in the variables j/j (typically exchanges or sign 
changes), but the equations of motion are invariant un- 
der these transformations. Thus, if a certain flow has 
this symmetry, it leads to constraints on the variables 
yi, and if it does not have this symmetry immediately a 
new flowfield can be obtained by applying this symmetry 
transformation. We do not attempt to analyze the full 
symmetry structure here and confine our discussion to 
two illustrative examples which are relevant for the sta- 
tionary states discussed below. Demanding invariance of 
the flow field to the reflection symmetry T\ leads to the 
following constraints on the variables yf. 



yi = V3 = yh = ys = 2/is = 
ye = vt 2/10 = -yu 

2/11 = 2/13 2/16 = -2/18 2/17 = 2/19- 



(25) 



The non vanishing components, 2/2 , 2/4) 2/6 = 2/7 1 2/9) 
2/10 = -2/12, 2/11 = 2/13, 2/14, 2/16 = -2/18, 2/17 = 2/19 thus 
define a 9 dimensional subspace. 

For the combined symmetry T1T2 we find the con- 
straints 



We achieved the impenetrability of the plates by re- 
quiring the mirror symmetry: 



(x,y,-z) 



u y (x,y,z) . 



(20) 



The reduction from 36 to 19 modes was achieved by 
restricting the dynamics to a subspace where the flow has 



the point symmetry around x 
the middle of the shear layer, 



(0,0,7r/2) 2 , a point in 



u y (x,y,z + ir/2) 



~ u v {-x,-y,-z + Tt/2) 



(21) 



In addition, there are further symmetries that can be 
used to reduce the phase space. There is a reflection on 
the y-z-plane, 



2/1 = 2/3 = 2/5 = 2/8 = 2/15 = 
2/6 = -2/7 2/10 = 2/12 

2/11 = -2/13 2/16 = 2/18 2/17 = -2/19 



(26) 



and again a 9 dimensional subspace with non vanishing 
components y 2 , 2/4, 2/6 = ~2/7, 2/9, 2/io = 2/12, Vn = -2/13, 
2/14, 2/16 = 2/18, 2/17 — —2/19- The dimensions of the invari- 
ant spaces vary from a minimum of 6 (for each a T\T$ 
and T1T2T3 invariance) and a maximum of 10 (for T2T3- 
invariance) . 

As mentioned, one can classify flows according to their 
symmetries. The most asymmetric flows are eightfold 
degenerate as the application of the eight combinations 
of the symmetries give eight distinct flows. The laminar 
flow profile is invariant under all the linear transforma- 
tions and is the only member of the class with highest 
symmetry. The other stationary states discussed below 
fall into equivalence classes with eight members or four 
members if they are invariant under T\ or T\T 2 . 
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III. DYNAMICS OF PERTURBATIONS 

A stability analysis shows that the laminar flow profile 
is linearly stable for all Reynolds numbers. The matrix 
of the linearization is non-normal with a block structure 
along the diagonal. To bring this structure out more 
clearly, it is best to order the equations in the sequence 
1, 2, 3, 4, 5, 7, 15, 8, 9, 14, 13, 19, 12, 18, 6, 11, 17, 10, 
16. The matrix of the linearization then is upper diago- 
nal, with a clear block structure: there are 10 eigenval- 
ues isolated on the diagonal, three 2x2 blocks and one 
3x3 block as well as several couplings between them in 
the upper right corner. While some eigenvalues can be 
complex, all of them have negative real part as shown in 
Fig. ||. For vanishing viscosity, the eigenvalues become 
zero or purely imaginary. 

Large amplitude perturbations, however, need not de- 
cay. Already in the linear regime the non-orthogonality 
of the eigenvectors can give rise to intermediate amplifi- 
cations into a regime where the nonlinear terms become 
important P-|l0|]. In a related study on plane Couette 
flow jl6) we used the lifetime of perturbations to get in- 
formation on the dynamics in a high-dimensional phase 
space. As in that case, the amplitude of the velocity field 
in the z-direction indicates the survival strength of a per- 
turbation. Linearizing the equations of motion around 
the base flow Uo gives for the perturbation u' the equa- 
tion 

d t VL = -(u • V)u' - (u' ■ V)u - Vp' + i^Au' . (27) 

The second term on the right hand side describes the en- 
ergy source for the perturbation, and depends, because 
of Uo = uq(z) e y and thus 



(u' ■ V)u = u' z d z u (z)e y 



(28) 



in an essential way on the ^-components of the pertur- 
bation. Thus, if the amplitudes y%, yio, y\%, yu, yie 
and 2/18 become too small, the decay of the perturbation 
cannot be stopped any more. These modes account also 
for most of the off-diagonal block-couplings. A model for 
sustainable shear flow turbulence has to include some of 
these modes. 

We chose a fixed initial flow field with a random selec- 
tion of amplitudes yx, ■ ■ ■ , yxa, scaled it by an amplitude 
parameter A and measured the lifetime as a function of 
A and Reynolds number Re. Fig. || shows the time evo- 
lution of such a perturbation at Re — 400 with one mode 
driven and for different amplitudes. For small A there is 
an essentially exponential decay, whereas for larger am- 
plitudes the perturbation swings up to large amplitude 
and shows no sign of a decline at all. The results for 
many amplitudes and Reynolds numbers are collected in 
Fig. H in a landscape plot. For small Reynolds number 
and/or small amplitude the lifetimes of perturbations are 
short, indicated by the light areas. For Reynolds num- 
bers around 100 isolated black spots appear, indicating 



the occurence of lifetimes larger than the integration time 
(which increases with Re so that t max /Re — 47r). The 
spottiness for Re between about 100 and 1000 is due to 
rapid changes in lifetimes from pixel to pixel. For Re 
above 1000 the long lifetimes dominate. These results 
are in good agreement with what has been observed in 
plane Couette flow. Fig. ||b shows a similar plot for the 
case with two modes driven; it is qualitatively similar, 
but quantitatively shifted to higher Reynolds numbers. 

In connection with the non-normality of the linearized 
eigenvalue problem it has been argued that the upper 
limit on the size of perturbations for which the non- linear 
terms in the dynamics can be neglected decreases alge- 
braically like Re~ a . Different exponents have been pro- 
posed, ranging from 1 to 3 || |ig , [L9| . It seems that for 
large Re (where the model is actually less reliable be- 
cause of the limited spatial resolution) the envelope of 
the long lived states in the fractal life time plot decays 
like Re- 1 . 

The sensitive dependence of lifetimes on initial condi- 
tions and parameters is further highlighted in Fig. ^ and 
|?]. The first shows the lifetime in the plane of the am- 
plitudes yi6 and yn at Reynolds number Re = 400 with 
all other components fixed. There is considerable struc- 
ture on many scales. One notes 'valleys' of short life- 
times between 'plateaus' of longer lifetimes and a granu- 
lar structure within both. The striations are reminiscent 
of features seen near fractal basin boundaries |20|. Fig. |?j 
shows successive magnifications of lifetime versus ampli- 
tude plots at Re=200. Even after a magnification by 10 7 
there is no indication of a continuous and smooth varia- 
tion of lifetime with amplitude. 



IV. STATIONARY STATES 

Motivated by the observation of new stationary struc- 
tures in plane Couette flow for sufficiently high Reynolds 
number [p~2|~[l~4| we searched for non-trivial stationary so- 
lutions and studied their generation, evolution and sym- 
metries. 

We computed the stationary states with the help of 
a Monte Carlo algorithm. The initial conditions for the 
j/i's were chosen randomly out of the interval [—1/2, 1/2] 
and the Reynolds number was chosen randomly ma- 
trixwith an exponential bias for small Re in the interval 
[10,10000]. With these initial conditions we entered a 
Newton algorithm. If the Newton algorithm converged, 
we followed the fixed point in Reynolds number as far 
as possible. We included about 200000 attempts in the 
Monte Carlo search. 

The stationary states found for a single driven mode 
are collected in Fig. |[ No stationary states (besides the 
laminar profile) were found for Reynolds numbers below 
about 190. Between 190 and about 500 there are eight 
stationary states which divide into two groups of four 
symmetry related states each. With increasing Reynolds 
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number more and more stationary states are found and 
they reach down to smaller and smaller amplitude. The 
envelope of all states reflects the i?e _1 behaviour found 
for the borderline where nonlinearity becomes important. 
For two driven modes (Fig. |^) the situation is similar. 

The appearance of the branches of the stationary states 
and in particular their coalescence near Re = 190 sug- 
gests that the states are born out of a saddle-node bi- 
furcation. And indeed, the eigenvalues as a function of 
Re show two eigenvalues moving closer together and col- 
lapsing at zero for Re = 190.41 (Fig. [To]). However, these 
eigenvalues are not the leading ones, so that one set of 
states has three unstable eigenvalues, the other two un- 
stable ones. It is thus a 'saddle-node' bifurcation into 
unstable states. 

With increasing Re more and more stationary states 
appear, partly through secondary bifurcations, partly 
through additional saddle-node bifurcations. Their num- 
ber increases rapidly with Reynolds number (Fig. |Tl| ) 
and this increase goes in parallel with the increase in 
density of long lived states, Fig. [| The detailed struc- 
ture of the bifurcation diagram is rather complex and 
has not yet been fully explored. We note here that the 
various stationary states may be grou ped according to 
their symmetries introduced in section IIP and that we 
found only stationary states which belong to equivalent 
classes with four or eight members. The stationary states 
of the classes with four members are invariant under the 
transformation T\ or T\T^. In addition, there are forward 
directed bifurcations generating two new branches with 
the same symmetry properties (eight or four member 
class) and inverse bifurcations of two branches belong- 
ing to eight member equivalent classes. We also found a 
backward directed bifurcation generating branches of an 
eight member class, which is born out of a four member 
class branch. The scenarios described above are marked 
in the bifurcation diagram Fig. g[ 



V. CONCLUDING REMARKS 

The few degrees of freedom shear model introduced 
here lies halfway between the simplest models of non- 
normality and full simulations. Its dynamics has turned 
out to be surprisingly rich. There are a multitude of 
bifurcations introducing new stationary states besides 
the laminar profile, there are secondary bifurcations, and 
the distribution of life times shows fractal structures on 
amazingly small scales. It seems that as one goes from 
the low-dimensional models via the present one to 
full simulations one notes not only an increase in numer- 
ical complexity but also the appearence of qualitatively 
new features pl| . 

The simplest models with very few degrees of free- 
dom focus on the non-normality of the linearized Navier- 
Stokes problem and emphasize the amplification of small 
perturbations. If the non-linearity is included a transi- 



tion to another kind of dynamics, sometimes as simple 
as relaxation to a stationary point, is found |l9[ . 

Next in complexity are models like the one presented 
here that share with the few degree of freedom models 
the amplification and the transition but the additional 
degrees of freedom allow for chaos. When nonlinearities 
become important the dynamics does not settle to a fixed 
point or a limit cycle but continues irregularly for an es- 
sentially unpredictable time. The time is unpredictable 
because of the fractal life time distribution which seems 
to persist down to amazingly small scale: tiny variations 
in Reynolds number or amplitudes of the perturbation 
can cause major variations in life times. This fractal be- 
haviour is the new quality introduced by the additional 
degrees of freedom. Indications for this behaviour are 
seen in the experiments by Mullin on pipe flow p^| . It 
is interesting to ask just how few degrees of freedom are 
necessary to obtain this behaviour. Reducing our model 
to the Ti subspace gives one with just 9 degrees of free- 
dom (comparable in number and flow behaviour to the 
ones of Walcffc [Q) that still shows this fractal life time 
distributions. Further reduction, as in the four mode 
model of po[ , seems to eliminate them. 

The full, spatially extended shear flows share essential 
features with the model but add new problems. Spatially 
resolved simulations of the present model JT||] as well as 
pla ne Couette flow with rigid-rigid boundary conditions 
jfi^fTit show the occurence of additional stationary states 
at sufficiently high Reynolds number that are unstable. 
A novel and as yet unexplained feature in spatially ex- 
tended plane Couette flow, which we believe to be con- 
nected to the high dimensionality of phase space, is the 
difference between Reynolds numbers where the first sta- 
tionary states are born (about 125 in units of half the 
gap width and half the velocity difference) and the ones 
where experiments begin to see long lived states (about 
300-350) ||. 

The fractal life time distributions have obvious simil- 
iarities to chaotic scattering |m| , [25|]2"c| 1 . Drawing on this 
analogy one would like to identify permanent structures 
in phase space away from the laminar profile that could 
sustain turbulence. This has partly been achieved by the 
search for stationary states. Many have been found but 
irritatingly only for Reynolds numbers above about 190 
while long lived states seem to appear much earlier. The 
solution to this puzzle must be periodic states and in- 
deed we have found a few periodic states in a symmetry 
reduced model at lower Reynolds numbers, close to the 
occurence of the first long lived states. This suggests that 
the dynamical system picture that long lived states have 
to be connected to persistent structures in phase space is 
tenable. 

There are several features of the model that can be 
studied further. In particular, quantitative characteriza- 
tions of the fractal life time distribution, visualizations 
of the flow field, a detailed analysis of the primary and 
secondary bifurcation, an investigation of the dependence 
on the aspect ratio of the periodicity cell are required and 



G 



look promising. We expect the lessons to be learned from 
this simple model to be useful in understanding the dy- 
namics of full plane Couette and other shear flows. Work 
along these directions continues. 
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FIG. 1. The 19 wave vectors ki . . . kig. The full set is ob- 
tained by complementing with — k 4 . Thus, there are only 
three vectors on the symmetry plane z — 0, six each on the 
two levels above and a single one on the third plane with 
k z = 3. 
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FIG. 2. The laminar profile in case of one (a) or two driven 
modes (b), compare eq. [n]. 



7 





0.0 0.2 0.4 0.6 0.8 1.0 



time [Re*d/u ] 

FIG. 4. The dynamics of perturbations for one driven mode 
at Re = 400. The perturbation was selected randomly and 
scale by factors 3, 5, 7, 9 and 11, from bottom to top. 




FIG. 5. Lifetime of perturbations as a function of ampli- 
tude and Reynolds number for the case of one driven mode 
(a) and two driven modes (b). The black regions correspond 
to lifetimes larger than T — 4n ■ Re, the white regions to life- 
times shorter than T/10 The grey levels interpolate linearly 
between these levels. 
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FIG. 6. Magnification of the fractal landscape of lifetimes 
as a function of the amplitudes yie and yi7 for the same per- 
turbation as in Fig. at Re = 400. 
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FIG. 7. Lifetimes of perturbations as a function of am- 
plitude for the case of one driven mode at Re = 200 and 
successive magnifications by a factor of 10. 
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FIG. 8. Stationary states for a single driven mode (a) and 
a magnification (b) near the leading saddle node bifurcation 
near Re =190. 
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FIG. 9. Stationary states for two driven modes. Compared FIG. 11. Proliferation of stationary states for one and two 

to Fig. ^ there seem to be more states and the next bifurcation driven modes, 
is a lot closer to the leading one. 
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FIG. 10. Eigenvalues of the two branches of the station- 
ary states a and b near the saddle-node bifurcation around 
Re ~ 190. Note that indeed two eigenvalues with real positive 
and negative real parts at zero, but there are also eigenvalues 
with positive real part. 
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